Center for Turbulence Research 
Proceedings of the Summer Program 1992 


307 


N94-rl^76S 

~ P 

DNS and modeling of the interaction 
between turbulent premixed flames and walls 

By T. J. Poinsot 1 AND D. C. Haworth 2 


The interaction between turbulent premixed flames and walls is studied using 
a two-dimensional full Navier-Stokes solver with simple chemistry. The effects of 
wall distance on the local and global flame structure are investigated. Quenching 
distances and maximum wall heat fluxes during quenching are computed in laminar 
cases and are found to be comparable to experimental and analytical results. For 
turbulent cases, it is shown that quenching distances and maximum heat fluxes 
remain of the same order as for laminar flames. Based on simulation results, a 
‘law-of-the-wall’ model is derived to describe the interaction between a turbulent 
premixed flame and a wall. This model is constructed to provide reasonable behavior 
of flame surface density near a wall under the assumption that flame-wall interaction 
takes place at scales smaller than the computational mesh. It can be implemented in 
conjunction with any of several recent flamelet models based on a modeled surface 
density equation, with no additional constraints on mesh size or time step. 


1. Introduction 

The understanding and modeling of turbulent phenomena that occur near walls 
is a formidable task. Even in the absence of chemical reaction, building ‘law-of- 
the-wall’ models or low-Reynolds-number models is an ongoing research subject 
and no satisfactory practical solution is yet available for general use in engineering 
codes. The situation is even more difficult when a flame is present. Combustion 
is strongly influenced by the presence of walls which may cause flame fronts to 
quench, for example. Moreover, the flame has a significant effect on the flow in 
the vicinity of the wall as well as on the heat flux to the wall. For these reasons, 
modeling of flame-wall interactions in turbulent situations is an important issue. 
Still, few experimental or modeling results have been reported, and most present 
models for turbulent premixed combustion do not use any specific corrections for 
near-wall effects. At best this may result in errors in the prediction of the reaction 
rate and of the wall heat fluxes and temperatures. In some cases, the absence of 
any reasonable approximation of the wall effects leads to numerical difficulties and 
to the use of ad-hoc numerical corrections to obtain solutions. These corrections 
are based on pragmatic rather than on physical grounds. Our objective here is to 
explore the flame-wall interaction mechanisms at a fundamental level using direct 
numerical simulations (DNS). The understanding thus obtained provides a sound 
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basis for a model which can be viewed as a law-of-the-wall approach for turbulent 
premixed combustion. 


2. Flame— wall interaction in laminar flows 

Before considering turbulent cases, it is useful to consider results obtained on wall 
quenching of laminar flames (Jarosinski 1986, Huang et al. 1986). For these flows, 
two important quantities have been introduced: the minimum distance between the 
flame and the wall 6q and the maximum heat flux 4 >q through the wall during 
the interaction with the flame. Most authors normalize the wall distance y by a 
characteristic flame thickness d = Aj/(piCps{*) and define the local Peclet number 
to be P e = yjd (Here a subscript ‘1’ refers to reference properties in the fresh 
gases). Therefore, the quenching distance is often expressed by its Peclet number 
PeQ- PtQ = S Q /d (Huang et al. 1986, Vosen et al 1984, Lu et al. 1990). The wall 
heat flux may be normalized by the laminar reference ‘flame power’ (heat release) 
to yield <j> = 4 , q/(pi YpS° AH) where pi and Yp designate the fresh-gas density and 
fuel mass fraction, s® is the unstretched laminar flame speed, and A H is the heat 
of reaction (YpAH = c p (T 2 — Tj) if T\ is the temperature of the fresh gases and 
T 2 is the adiabatic flame temperature). These two quantities may be correlated in 
laminar flows. If one assumes that at quenching, the wall heat flux 4 >q is due to 
heat conduction in the gas layer of thickness 6q, one can write 




A 


T 2 -T w 

6q 


( 1 ) 


where A is the gas thermal conductivity and T w is the wall temperature. From 
Eq. (1), we obtain a relation which is expected to hold when the wall heat transfer 
is dominated by diffusion: 


_ T 2 - Tu, 

T 2 - T\ P c q 


Three typical situations have been studied in the past (Figure 1): 

1. Head-on quenching 

When a flame front reaches a cold wall (T w = T x ) at a normal angle, head-on 
quenching (HOQ) occurs (Figure la). This case has been studied numerically and 
experimentally (Huang et al. 1986, Jarosinski 1986, Vosen et al. 1984). Results 
suggest that quenching occurs for Peclet numbers of the order of three. Heat flux 
measurements indicate values of <j> of the order of 0.34, which is consistent with the 
value predicted by Eq. (2) with P e Q ss 3. In terms of simple physics, this result 
that a flame stops propagating towards the wall when the heat losses to 
the wall are equal to about one-third of the nominal flame power. The fact that (j> 
is almost constant for different fuels (Huang et al. 1986) suggests that the problem 
is thermally controlled and that simple chemistry may be used to compute this 
phenomenon. 
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Figure 1. 


Configurations for flame-wall interaction studies in laminar flows. 


2. Side-wall quenching 

When a flame propagates parallel to a wall, the situation is different. Only localize 
quenching of the flame near the wall occurs (Figure lb). T ^ ls sltuatlon has ^ en 
studied theoretically by Von Karman & Millan (1953) and Makhvi adze & Mehkov 
(1991), and experimentally by Lu et al. (1990) and Clendemng et al. (1981). Peclet 
numbers in this case are of the order of seven suggesting values for <t> of about 0.16 
fEq (2)) Asymptotic theories of non-adiabatic flames also may be used to pre ic 
the quenching distance (Williams 1985): these predict the same order of magnitude 

for P e Q. 

3. Tube quenching 

Total flame quenching may occur in a tube if its diameter is sufficiently small 
(Lewis k Von Elbe 1987, Jarosinski 1986, Fairchild et al. 1984) (Figure lc). This 
phenomenon is exploited, for example, in the design of flame arrestors: these are 
ensembles of tubes with diameters smaller than the quenching distance so tha a 
flame cannot propagate through them. Peclet numbers in this < «*e are ' based on 
the tube diameter and are of the order of 50 (Aly k Herman ce 1981). We will not 
consider this configuration here since in most practical situations, the dimensions 
of the system (e.g., the size of the combustor chamber) are too large to induce total 

quenching. 
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3. Numerical method and configuration 

To study numerically how a laminar or a turbulent flame interacts with a cold wall 
we have utilized a DNS code developed at Stanford (Poinsot et al. 1991, Poinsot 
&: Lele 1992). Since this code has been described elsewhere, we will recall only 
its principal features here. We consider a compressible viscous reacting flow. The 
chemical reaction is represented by a single-step mechanism, R (reactants) -*■ 
P (products) where the reaction rate w R is expressed as, 


wr = BpY R exp 



( 3 ) 


This can be interpreted as a binary reaction where one of the reactants (Y R ) is 
always deficient. Following Williams (1985) we cast this expression in the form, 


»* = A,y R exp( r ^L-| 5 ). (4) 

Here 0 is the reduced temperature 0 = (T - 7) )/(T 2 - T,), and T 2 is the adiabatic 
flame temperature. The activation temperature is T a and the coefficients A, a, and 
P are, respectively, the reduced pre-exponential factor, the temperature factor,' and 
the reduced activation energy, 


A = B exp(-(3/a), a = (T 2 - T x )/T 2 , and /? = aT a /T 2 . (5) 

Fluid properties follow the equations, 

P = Pi{pTi/ P iT), p = p i (T/T 1 ) i , 

Le = \j pDcp = constant , Pr = p Cp /X = constant , (6) 

where p, A, and D are molecular diffusivities of momentum, internal energy, and 
species, respectively, and b is a constant. Using these assumptions and a Cartesian 
frame of reference, the conservation equations for compressible flows are solved using 
a high-order finite difference scheme (Lele 1992). 

The calculations are initialized with reactants on one side of the computational 
domain and products on the other; these are separated by a laminar premixed 
flame. The wall is located on the reactant side of the domain (Figure la). All 
velocity components are zero on the wall and the wall temperature is imposed. For 
all the cases shown here, the wall temperature is equal to the fresh gas temperature, 
T\- On lateral boundaries, periodic conditions are enforced. On the right-hand side 
of the domain, non— reflecting boundary conditions are used (Poinsot Si Lele 1992). 

The initial velocity field (turbulence spectrum) and spatial distribution of reac- 
tant mass fraction are specified at t ~ 0: the system is then allowed to evolve in 
time. The initially planar flame is convected and strained by the turbulence while 
the combustion influences the fluid mechanics through dilatation and temperature- 
dependent properties (Eq. (6)). After some time (typically 5 to 20 flame times in 
these computations), the flame reaches the wall and begins to interact with it. 
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4. Diagnostics 

In the present case, we are especially interested in the effect of the wall on the 
flame structure. Postprocessing of the two-dimensional computed fields (snapshots 
at fixed times) begins by defining a flame front as an isocontour of temperature 
T Once the flame front has been located, the local normal and local flame cur- 
vature are readily computed; curvature is taken to be positive for flame elements 
that are concave towards the products and conversely for elements concave towards 
reactants. One-dimensional cuts normal to the flame are taken: it is these pro- 
files that define the local ‘structure’ of the turbulent flame. We compare the local 
turbulent flame profiles with the steady one-dimensional laminar flame profile for 
the same chemistry and fluid properties. Of particular interest is the distribution 
along the flame of the normalized local flame speed (‘flamelet speed ) s n defined 
by 3n _ J dn / s°, , that is, the integral of the reaction rate profile in a direction 

locally normal to the flame. 

The fixed chemical parameters used for this study are summarized m Table 1. 1 1 e 
flame speed and thickness are normalized respectively by the sound speed c and by 
the reference length d = A 1 /(p 1 c p s?). The flame thermal thickness 6, is based on 
the maximum temperature gradient: 6 ? = (T 2 - Ti)/Max( dn ). 


Table I. Fixed parameters for DNS of flame-wall interaction. 


a 

0 

A 

b 

Pr 

Le 

s°/c 

6°/d 

0.75 

8.00 

146. 

0.76 

0.75 

1.0 

0.016 

3.8 


5. DNS of the interaction between a laminar flame and a wall 

To check the accuracy of the model, laminar runs were performed first Figure 
2 presents time variations of the flame distance to the wall as well as the flame 
power Pl s n c p (T 2 - T,) and the normalized wall heat flux <f>. There time has been 
normalized by the flame time t F = 6°/s° and y is the distance from the wall. The 
values obtained from DNS for this case are P eQ = 3.4 and <j> - 0.36. These values 
are in good agreement with experimental data (Lu et al. 1990, Vosen et al. 1984) 
and with the simple model given by Eq. (2). Although total quenching occurs at 
Peo = 3-4, the flame senses the presence of the wall before this time: the flame 
speed s„ begins to decrease when the wall distance is less than a distance S T given 
by p t _ s T /d sa 8. Therefore, two zones are necessary to describe the near-wall 

*(i) The ‘quenching’ zone stretches from the wall to a local Peclet number y/d of 
about 3.5 (0 < y < 6 q). In this zone, no reaction ever takes place. 

(ii) The ‘influence’ zone goes from the wall to a Peclet number y/d of about 8 
(0 < y < S T )- Any flame entering the influence region will start sensing the wall 
and will eventually get quenched. The time t Q needed for the flame to quench after 
it enters the influence zone is of the order of two flame times, t Q = 2 t F - 26, /s, . 
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FIGURE 2. Numerical results for a laminar flame- wall interaction. 

6. DNS of the interaction between a turbulent flame and a wall 

The interaction between the wall and a turbulent flame front is characterized by 
three effects: 

(i) A local thermal effect, by which heat losses to the wall affect the flame structure 
and result in local quenching. 

(ii) A geometrical effect which limits the spatial extent of the flame brush and 
reduces the flame-brush size in the vicinity of the wall. 

(iii) A laminarization effect which is an indirect effect of the wall on the flame. The 
wall affects the structure of the turbulence and leads to laminarization immediately 
adjacent to the wall. This induces a strong decrease of the turbulent stretch and 
thereby a decrease of the flame area. 

Preferential species diffusion is beyond the scope of the present simple-chemistry 
investigation. 

The structure of turbulence near the wall is clearly an important issue in the 
latter two questions. The configuration studied here corresponds to the ‘shear-free’ 
boundary layer in which turbulence with no mean shear interacts with a wall. This 
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may be realized experimentally in a wind tunnel, for example, using walls that move 
at the mean flow speed. The shear-free boundary layer has been studied experi- 
mentally by Thomas k Hancock (1977) and Uzkan k Reynolds (1967), th^retically 
by Hunt k Graham (1978), and numerically by Biringen k Reynolds (1981). I he 
flow structure can be summarized as follows. Starting from initially isotropic homo- 
geneous turbulence, the wall induces a perturbation zone whose thickness increases 
with time. In this zone, viscous effects are important, and the turbulence is damped. 
Moreover, by imposing zero normal velocities, the wall increases velocity pertur a- 
tions in planes parallel to the wall (at sufficiently high Reynolds numbers). The 
role of this near-wall turbulence structure on flame evolution is difficult to quantify 
in the present simulations. We will concentrate instead on the thermal effects. 

The parameters for the simulations reported here are summarized in Table 11. 
There u' is the rms turbulence velocity, L t is the length scale of the most energetic 
vortices in the initial turbulence spectrum, and / is the turbulence integral scale 
based on two-point velocity correlations. The initial turbulence field was chosen 
to produce small-scale turbulence near the wall and to impose zero velocity fluc- 
tuations at the wall. More sophisticated (three-dimensional) approaches will be 
necessary to produce more general results. 

Table II. Initial conditions for DNS of turbulent flame-wall interaction. 


Case 

“7*1 

Li/S] 


Rei = u'l/u 

II 

*4 

2D3 

6.25 

8.9 

2.85 

90 

280 

2D4 

6.25 

4.5 

1.43 

45 

140 

2D6 

6.25 

1.9 

0.64 

19 

60 


Figures 3 and 4 present snapshots at one instant in time of reaction rate and 
vorticity fields during the interaction between the turbulent flame and the wall. 
Time has been normalized by the flame time t F = $“/«“, and the maximurn value 
of the vorticity modulus is normalized by the characteristic flame strain s, /6, . The 
structure of the vorticity field is affected both by the flame (viscosity in the burnt 
gases dissipates vorticity rapidly) and by the wall (the normal velocity component 
close to the wall goes to zero while the parallel component increases). Vortex pairs 
appear to play a dominant role. In Figures 3 and 4, for example, one vortex pair 
attracts a part of the flame front towards the wall while another pair pushes a 
different part of the flame away from the wall (sequence t/t F = 4.4 to t/t F = 6-6)- 
This is confirmed by Figure 5 which presents the time evolution of the minimum 
and maximum wall distances in terms of Peclet numbers. When the first flame 
element touches the wall at time t/t F ~ 6, the most distant element is moving 
away from the flame front. This ejection of flame elements away from the flame 
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Vorticity Reaction rate 




FIGURE 3. Vorticity and reaction rate contours for a turbulent flame- wall in- 
teraction (Case 2D4). Solid lines denote clockwise vorticity, dashed lines denote 
counterclockwise. 


front may induce local counter-gradient turbulent diffusion of flame surface density 
which might be important in modeling. 

As far as the thermal effect is concerned, results obtained during this simulation 
and during other simulations of the same type lead to a simple result (Figure 6): the 
maximum local heat flux to the wall corresponds within 10 percent to the laminar 
heat flux (<f> = 0.36), and the quenching distance 6q is equal to the value obtained in 
laminar cases (P £ q = 3.4). Although the initial conditions used for these simulations 
lead to large velocity perturbations near the wall, the heat flux to the wall appears 
to be controlled mainly by heat diffusion, and the local instantaneous maximum 
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Vorticity 






FIGURE 4. Vorticity and reaction rate contours for a turbulent flame wall inter 
action (Case 2D4, continued). 


heat flux obtained during quenching is the same as for the laminar cases. 

Figure 7 presents a scatter plot of normalized flamelet speed s n versus distance 
to the wall. This plot exhibits different behaviors for different flame elements 
(‘flamelets’) as they approach the wall. No flame elements approach closer than 
the laminar quenching distance 6q to the wall, and the minimum Pec let number 
is equal to about 3.5. Branch 1 corresponds to flames reaching the wall at normal 
angles (head-on quenching): these follow the curve (solid line) predicted by the lam- 
inar HOQ computation (Section 5), quenching at a Peclet of close to 3.5. Branch 
2 corresponds to flamelets which disappear at Peclet numbers of about 7. Exam- 
ination of DNS fields for these points suggest that these flamelets are propagating 
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Figure 5. Maximum and minimum (over all computational cells adjacent to the 
wall) flame-wall distances for Case 2D4. 



Time / Flame time 


Figure 6. Minimum, maximum, and mean (over all computational cells adjacent 
to the wall) wall heat fluxes for Case 2D4. 


parallel to the wall and not towards the wall, thus corresponding to the side— wall 
quenching situations described in Section 2. Branch 3 corresponds to flamelets 
which burn faster (accelerate) as they approach the wall, but subsequently quench 
on reaching a Peclet of about 3.5. At this point, no explanation for Branch 3 is 
proposed. The number of flamelets following Branch 2 is small compared to the 
other branches: most flamelets reach the wall at a normal angle. This is confirmed 
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FIGURE 7. Flamelet consumption rate s n versus wall distance at time t ft p 4 
(Case 2D6). 

by the distributions of flame curvature presented in Figure 8. While the flame is 
still far from the wall (t/t F = 0.5), the mean curvature conditioned on distance 
to the wall is almost symmetrical about zero. As the flame brush approaches the 
wall (t/tF = 4). flamelets flatten and curvature diminishes. Furthermore, positive 
values of curvature are clipped more strongly than negative values indicating that 
flamelets are predominantly concave towards the fresh gases and that most of them 
will reach the wall at close to a normal incidence angle. 

The distance at which flamelets begin to sense the presence of the wall (i.e., 
where their local flamelet speed begins to drop) is given by a Peclet number of the 
order of 10 (Figure 7), close to that found for the laminar simulations. Although 
the existence of the three different branches suggests a more complex pattern than 
simple head-on quenching, it appears that an influence zone may be defined for 
turbulent cases whose thickness (in Peclet number) is something close to 10. 

It appears that both the quenching zone thickness (Sq) and the influence zone 
thicknesses ( S T ) have similar values for turbulent and laminar premixed flames. This 
has some important consequences. Consider, for example, a reacting boun ary 
layer. Invoking the usual normalizations, we denote by a superscript a wall- 
units-scaled value: y + = yu T /v, where u T is the friction velocity and i/ is the 
kinematic viscosity. The edge of the quenching zone is located at S Q - u T bq/v - 
tp 0 / Pr)(u T / s?). For most practical situations, u r (1 to 5 m/s) is of the order of 
the flame speed s° (0.3 to 1.5 m/s) so that 6+ is of the order of 1 to 10. That means 
that the quenching zone is located inside the viscous layer. Flamelets travelling 
from the free stream towards the wall first encounter laminar flow m the viscous 



318 


T. J. Point ot & D. C. Haworth 


region and only later are quenched. This result will be used in the modeling of the 
flame- wall interaction (Section 7). 



FIGURE 8. Mean curvature (conditioned on wall distance) for a flame far from the 
wall {t/tf = 0.5) and a flame reaching the wall {t/tp = 4) (Case 2D6). 

7. A law-o f—t he— wall for turbulent premixed combustion 

From the previous results, it is possible to construct law— of— the— wall models for 
premixed turbulent flames which we will designate here as ‘FIST’ models (Flame 
Interacting with Surface and Turbulence). This model will be derived in the frame- 
work of flamelet models: the dependent variable which will be modeled is the flame 
surface density (surface-to-volume ratio) £ as defined by Pope (1988). Implemen- 
tation is discussed in the context of a finite-volume method, although the concept 
may be applied to other numerical approaches. 

We consider a generic flamelet model in which the flame surface density £ evolves 
according to, 


5£ dU ] •£ 

dt dxi 


d?i 

-^ + S-D-Dq. 


( 7 ) 


This equation includes transport, turbulent diffusion {T t ), and source (5) and con- 
sumption terms ( D , Dq). With the exception of the Dq term which represents 
thermal quenching due to the wall, several recent flamelet models can be cast into 
this form (Cant et al. 1990, Candel et al. 1988, Boudier 1992, Cheng & Diringer 
1991). The development of the present FIST model is independent of the exact 
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Figure 9. Principle of the FIST model (law-nf-the-wall model for premixed 
turbulent combustion). 

form of the 5 and D terms. The assumptions which will be used are the following 
(see Figure 9): 

- We will consider a situation where a law-of-the-wall approach is used to de- 
scribe near-wall turbulence since this is the case for most practical codes based on 
Reynolds-averaged mean equations. 

- The computational cell size adjacent to the wall Ay is larger than the quenching 
zone and larger than the zone over which the wall modifies the free-stream turbu- 
lence structure (typically, in the case of a turbulent boundary layer, the first grid 

point is located at a y + larger than 200). 

- Outside the quenching zone (y > 6t), no thermal quenching occurs as shown 

by DNS (see Section 6) so that Dq = 0. 

- Inside the quenching zone (0 < y < 6r)» we need to estimate the characteristic 
time t Q at which flamelets are quenched. For laminar cases, t Q is of the order of 
two flame times as evidenced by DNS for laminar flames (see Section 5). Despite 
the existence of the three branches shown in Section 6, we will assume that all 
flamelets entering the influence zone in the turbulent case are quenched on a time 
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scale tq — 2 tp. Therefore, Dq may be written as Dq = Eq/*q where is 
the mean flame surface density in the quenching layer. This thermal quenching 
effect is supposed to be stronger than the usual consumption terms (due to mutual 
annihilation) in the S equation so that D may be set to zero. Since there is no 
turbulence inside the quenching layer (see Section 6) we will also assume that there 
are no source terms in this zone (5 = 0). 

- Inside an ‘inhibition zone’ which extends from y = 0 to we assume that 

turbulence is affected strongly enough by the wall to reduce the flame stretch to zero. 
This assumption is not based on the present DNS but rather on the observation that 
turbulence must be strongly damped in this zone, therefore reducing the turbulent 
stretch. 

- In the rest of the first computational cell (yt < y + < Ay + ), the normal form 
of the flamelet model is used. 

- Finally, for the sake of clarity, we will consider a situation where the flow is 
homogeneous in planes parallel to the wall (i.e., S is a function of y only), and we 
will also neglect mean velocities normal to the wall. These last assumptions allow 
mean convective terms ( ) to be neglected in the finite-volume expression of 
Eq. (7). This assumption is not restrictive, and is invoked only for convenience in 
writing the resulting modeled equations. 

Under the previous assumptions, the equations of a FIST model may be derived 
by integrating Eq. (7) in two zones: the influence zone and the rest of the first 
computational cell. The mean flame surface densities in the influence zone and in 
the first computational cell, respectively, will be defined by: Eq = f* T £ dy , and 

=•*= 3^1 

In the quenching zone, terms 5 and D are small compared to the thermal quench- 
ing effect. When Eq. (7) is integrated between y = 0 and y = St with this assump- 
tion, the following conservation equation for Eq is obtained: 

3Eq F(y = St) Sq , fi . 

di 6 t t Q ' W 

Integrating Eq. (7) between y = 6 t and y = Ay provides an evolution equation for 
the average flame surface density E] in this zone: 

°§r = A^h r my - ^ y = - o . (9) 

Equations (8) and (9) form a closed set which provides the flame surface density 
in the quenching region (Eq) and in the first computational cell (Ei). These two 
equations state that there is a sink mechanism for flame surface in the first com- 
putational cell: flamelets diffuse towards the quenching zone (^F(y = St) term in 
Eq. (9)) and later get quenched in this zone on the time scale <q. 

8. An equilibrium formulation for FIST models 

Although Eqs. (8) and (9) may be solved under this form in finite-volume codes, 
it is interesting to propose a simpler model ‘Equilibrium FIST’ in which the flame 
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surface density in the quenching region Eq may be eliminated. In the Equilibrium 
FIST model, four additional assumptions are invoked: 

- The turbulent diffusion term Ti is written as JF; = where u t is a turbulent 

diffusion coefficient and S c is a turbulent Schmidt number. 

- The flame surface-density profile inside the quenching zone is supposed to 

exhibit strong spatial variations compared to the profile outside this zone.^Then 
the diffusion term F y (y = S T ) may be estimated by T y (y = S T ) - EQ t T 
where E is a model constant of order unity. . _ , 

- The quenching zone is assumed to be in equilibrium, i.e., diffusion balances 
dissipation in Eq. (8). This allows us to derive an explicit expression for the flame 
surface density in the quenched region as a function of the flame surface density in 
the first cell Ej: 

E Q = El- aQ - where a Q = E^TT • ( 10 ) 

+1 &C 0 T 

The parameter a q is proportional to the turbulent diffusivity v t normalized by the 
quenching time and distance.! 

- The size of the first cell is supposed to be sufficiently large compared to the 
quenching distance (Ay >> S T ) to neglect S T in the RHS of Eq. (9). 

Under these assumptions, the Equilibrium FIST model provides the following 
conservation equation for the flame surface density Ei near a wall: 


as, 

dt 


1 ^-/ an ^ t a Q i c/i \ __ n 


( 11 ). 


Wall corrections appear here only as an additional diffusion term towards the wall 
(second term on right-hand side of Eq. (11)) and as a correction of the turbulent 
stretch 5(1 — y+ /Ay). All other terms may be estimated by classical finite-volume 
methods. In this model, some constants may be set directly from the present DNS 
results: the influence distance is given by S T = P*d = where the Peclet 

number should be of order 10 (Sections 5 and 6), and the quenching time scale t Q is 
given by t Q = 2 t F = 2l° F /s^ (Section 5). The parameter y? n has not been estimated 
from the present DNS results but might be determined by using a three-dimensional 
boundary-layer code. Reasonable estimates for this quantity are of the order of 50. 
The turbulent diffusivity u t appearing in the above formula may be estimated using 
standard expressions for this quantity near walls. Further improvements of the 
model may be based on an expression for v t which would take into account the 
counter-gradient diffusion of E mentioned in Section 6. 

The simplicity of this formulation allows it to be used in conventional Reynolds- 
averaged multidimensional flow codes without additional constraints on time step 


! By using an eddy-viscosity concept near the wall (l^i — A u'S'j'), aQ may be interpreted as 
a ratio between turbulence velocity near the wall and a characteristic quenching velocity ( a Q = 

A s:s T /tQ h 
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or grid size. Although the approach simplifies the actual physics of flame-wall 
interaction, it represents a significant improvement over approaches which fail to 
account explicitly for the influence of the wall on the turbulent flame. The Equi- 
librium FIST model accounts for turbulent diffusion of flamelets towards the wall 
and quenching on a time scale which is given by DNS. It also accounts in a crude 
way for the laminarization effect of the wall up to a distance given by yf n . 

9, Conclusions 

Calculations of premixed laminar and turbulent flames interacting with isother- 
mal walls have been reported. Quantitative results have been presented illustrating 
the influence of distance from the wall on the local and global flame structure. 
For laminar cases, the computed minimum distance between wall and flame (the 
‘quenching distance’) and the maximum wall heat flux during quenching have been 
found to be comparable with available experimental and analytical results. For tur- 
bulent cases, it has been shown that quenching distances and maximum heat flux 
remain of the same order as for laminar flames. Correlations between wall distance 
and flame structure suggest that thermal effects are important only very close to the 
wall and that the wall acts as a strong sink term for flame surface density. Based on 
these DNS results, a model has been proposed to take into account the interaction 
between the turbulent flame and the wall. The equilibrium version of this model 
may be implemented in conventional finite-volume codes together with flamelet 
models based on modeled surface density equations. Further tests are necessary to 
assess its performance. 

An important extension of the FIST model would be the development of a model 
for wall heat flux. Such a model could be based on the knowledge of the flame 
surface density in the quenching zone and on the correlations between wall heat 
flux and flame position. This issue will be addressed in future work. 
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